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Abstract 

Regge's method for regularizing euclidean quantum gravity is applied to two 
dimensional gravity. We use two different strategies to simulate the Regge path in- 
tegral at a fixed value of the total area: A standard Metropolis simulation combined 
with a histogramming method and a direct simulation using a Hybrid Monte Carlo 
algorithm. Using topologies with genus zero and two and a scale invariant integra- 
tion measure, we show that the Regge method does not reproduce the value of the 
string susceptibility of the continuum model. We show that the string susceptibility 
depends strongly on the choice of the measure in the path integral. We argue that the 
failure of the Regge approach is due to spurious contributions of reparametrization 
degrees of freedom to the path integral. 
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1 Introduction 

It is still an unsolved problem how a quantum theory of gravitation can be formulated 
using an euclidean path integral. Well-known problems are the unboundedness of 
the euclidean action, the perturbative non-renormalizability of the Einstein-Hilbert 
action and finding a consistent regularization, for instance by discretizing space- 
time, such that the path integral can be studied non-perturbatively with numerical 
methods. Here we shall focus on the problem of finding a regularization of the 
euclidean path integral. 

Most studies so far have either used dynamical triangulation [|l], |^ or Regge 
calculus 0. For recent review articles we refer the reader to refs. @, |^. For an 
approach along a different line, see refs. 0]. In both cases the path integral is 
replaced by a summation over simplicial manifolds. These simplicial manifolds are 
built up from d dimensional simplices, triangles in d = 2 dimensions, which are glued 
together along their d — 1 dimensional subsimplices, the edges in two dimensions. In 
the dynamical triangulation approach one keeps the edge lengths fixed and replaces 
the integration in the path integral by a summation over simplicial structures with 
all possible connectivities and numbers of simplices. In the Regge method the edge 
lengths are the dynamical degrees of freedom, while the connectivities are usually 
kept fixed. Here the path integral is replaced by an integration over the edge lengths. 
Both formulations have appealing features, but are also subject to criticism. 

In both approaches the metric tensor field, which is the integration variable in 
the continuum path integral, is replaced by a simplicial structure. With dynamical 
triangulation, the Einstein-Hilbert action is automatically bounded from below; in 
the Regge approach one usually adds a curvature square term to make the action 
bounded 0. It is far from clear, however, how the continuum integration measure, 
which includes gauge fixing and Faddeev-Popov terms, has to be translated to a 
summation over these simplicial structures. 

With dynamical triangulation a natural choice is to simply sum over all triangu- 
lated surfaces, only weighted by the action. Since it is impossible to use dynamical 
triangulation to study small fluctuations around a flat space in perturbation the- 
ory, the only way to test this conjecture is through numerical simulations. This 
should reveal the existence of a second order phase transition at which renormalized 
quantities, such as the volume and curvature fluctuations, become independent of 
the details of the underlying discrete structure. Studies of the phase structure have 
given some evidence for the existence of such a critical point in four dimensions 
[^, Hy], but the situation is far from clear Jll], [1^. It has furthermore been ques- 



tioned if the Monte Carlo methods which are presently in use are ergodic |13[- The 
strongest support for the dynamical triangulation approach comes from its success 
in two dimensional gravity. Here it has been shown in detail, that the results of the 



continuum theory are reproduced [^, |, 0] 
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In the Regge approach it is less easy to choose a natural measure for the integra- 
tion over the link lengths. In contrast to the dynamical triangulation approach, it is 
possible to study the Regge path integral in perturbation theory. Here it has been 
shown that the path integral with a simple (scale invariant) measure, corresponds 
to the continuum path integral without gauge fixing or Faddeev-Popov terms |T^ . 
This suggests that also outside perturbation theory, the integration includes degrees 
of freedom which correspond to coordinate transformations; we shall discuss this in 
more detail below. To remove these spurious degrees of freedom, one would expect 
that the measure should include gauge fixing and Faddeev-Popov terms. One has 
often commented on this issue, but in actual simulations only the scale invariant 
measure or simple variations thereof have been used. 

Even though the justification of the measure used in dynamical triangulation 
appears to be equally obscure as in the Regge approach, this latter approach has 
been lacking the compelling success of the dynamical triangulation method in two 
dimensional gravity. This has motivated some time ago an investigation of the 
Regge approach in two dimensional gravity . This numerical study claimed that 
also the Regge method reproduces correctly the scaling of the partition function, in 
particular the so-called string susceptibility of two dimensional continuum gravity 
with central charge zero was found. The string susceptibility was not found to 
coincide with the continuum prediction for non-zero central charge. In the present 
work we describe the results of a similar study of the central charge zero model. 



which we believe improves several of the shortcomings of ref. [jT6[. We shall use 
the same Regge action, which includes an term, and the same scale invariant 
integration measure as in ref. ||16[. Unfortunately, we find wrong values for the 
string susceptibility, when the simplicial structure has the topology of a sphere or a 
"bi-torus" (a manifold with two handles). We will show explicitly that the scaling of 
the Regge partition function depends very sensitively on the details of the integration 
measure. 

The paper is organized as follows. In sect. 2 we review the continuum scaling 
relations for the path integral over surfaces with a fixed area. We discuss in some 
detail how the Regge skeletons can be viewed as discretizations of the metric field 
of continuum manifolds. Sect. 3 deals with the derivation of scaling relations for 
the Regge path integral over skeletons with a fixed area. We shall show that these 
scaling relations reveal a sensitive dependence of the string susceptibility on the 
integration measure. In sect. 4 we discuss the technical aspects of our simulation 
methods. As a novelty, a Hybrid Monte Carlo algorithm is introduced which can 
be used to simulate the fixed area Regge path integral. Our numerical results are 
presented in sect. 5 and a brief summary and some final remarks are contained in 
sect. 6. 
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2 Regge approach to two dimensional gravity 
2.1 Continuum scaling relations 

Let us consider the euclidean path integral for two dimensional pure gravity (the 
central charge zero model), 

S{g) = J d^x ,/g {X + kR + (3R^) . (2.2) 

In two dimensions the curvature term / d?x^R = Anx is a topological invariant 
(Gauss-Bonnet theorem) depending only on the Euler characteristic x = 2(1 — h), 
where h denotes the number of handles or genus of the surface. The genus is zero 
for the sphere, one for the torus and two for the bi-torus. We have included a 
term in the action with a coupling constant jS which has a dimension of 1/mass^. 
The coupling constant k, is dimensionless and the cosmological constant A has a 
dimension of mass^. The delta function in ( |2.1| ) imposes the constraint that the 
area of the surface is equal to A. Formally the measure is normalized by the infinite 
volume of the diffeomorphism group, which should cancel the infinite factor that 
arises from the integration over the gauge (reparametrization) degrees of freedom in 
the path integral. The continuum measure Vg is ill-defined because it involves an 
infinite product over all space-time points. 

On a formal level, i.e. ignoring the measure, one would expect that 

Z{A) oc A-^exp[-\A] , (2.3) 

if we set (3 in the action (|2.2|) equal to zero. However, after regularization and 
fixing the gauge using the usual Faddeev-Popov procedure, one finds that the mea- 
sure breaks the conformal (scale) symmetry and the A dependence of the partition 
function is given by 

Z{A) oc A'^^"^-' exp[-XRA] , (2.4) 

where Xr is the renormalized cosmological constant. The string susceptibility 7str 
depends on the genus h of the surface, 

7str = 2-^{l-h). (2.5) 



The relations (p.4|) and ( p.5|) were first derived for the special case h = using the 
light cone gauge and were later confirmed and generalized to topologies of 
arbitrary genus h using the conformal gauge. Quantum fluctuations are seen 
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to lead to a deviation from the naive scaling dimension of the partition function for 
h ^ 1 and renormalize the cosmological constant. 

For small, but non-zero (3 one expects that the result ( p.4D remains valid in 
the regime where A ^ because the term should only modify the short 

distance behavior of the theory. The R^ model for A//? — > has been studied 
in ref. [^|. Using the dynamical triangulation method it has been demonstrated 



numerically that the addition of a small R^ term did not affect the value of the 



string susceptibility M 21 



In order to study this two dimensional gravity model using the Regge regulariza- 
tion, one first has to express the action (|2.2| ) in terms of the edge lengths and then 
integrate over all length variables with an appropriate choice for the measure. In 
order to expose more clearly the role of reparametrization invariance we shall recall 
briefly that the Regge action can be viewed as a discretization of the continuum 
Einstein-Hilbert action. This establishes a link between the edge lengths of a Regge 
skeleton and the metric field in the continuum and furthermore sheds light on the 
problem of choosing the integration measure in the Regge approach. 



2.2 Regge discretization 

Usually the Regge skeletons are considered as special, piecewise fiat, continuum 
manifolds. Here we want to take a somewhat different point of view, in which a 
Regge skeleton is a discretized version of an underlying continuum manifold which 
is defined by a metric tensor field g^j_y{x). For simplicity we shall assume that the 
continuum manifold T has the topology of a torus. Then we can choose global 
coordinates {x = (xi, 0:2)} on this surface with metric g^y{x). The coordinates are a 
(globally defined) map from a fiat unit torus To = {x\xi G [0, 1), X2 G [0, 1)}, into the 
curved surface of T. We can construct a grid on the curved torus T by first defining 
the grid A = = n^ao; = 1, • • ■ , A^; oq = on To, as shown in fig. |I] and 

then mapping it onto the curved surface. In this way we can obtain triangulations of 
T where the triangles are spanned by the triplets of coordinates {(nia, n2a), ((ni + 
l)a,n2a), {{rii + l)a, (n2 + l)a)} and {(nio, n2a), (nio, (n2 + l)a), ((^i + l)a, (n2 + 
l)a)}, cf. fig. |I[ The edge lengths to be used in the Regge skeleton, can be defined 
in various ways: one can take the proper length of images of the corresponding edges 
on the unit torus Tq, one can use the geodesic distance or one can use the distances 
in some fiat embedding space. 

Using such a construction, one can give the relation between the continuum 
metric g^^ and the hnk lengths of the Regge skeleton. 



hi = y/gu do + 0{al), 
hi = v^"o + 0(ao). 
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Figure 1: Triangulation of the torus with N = A. The boundaries with the same 
labels have to be identified. The lattice spacing is denoted by a^. Each vertex in the 
skeleton has coordination number six. 

hi = \Jgii + 25fi2 + 5-22 flo + 0{al), (2.6) 

where = 1/N is the lattice distance on the unit torus Tq. The metric tensor 
g^^i, in ( p?6|) is evaluated at the site Xj = (niao,n2ao) G A and the three links 
lii, hi and hi are the edges on T that are associated with the coordinate pairs 
{{niao,n2ao), + l)ao, n2ao)}, {(niOo, n2ao), (niOo, (n2 + l)ao)} and {(niOo, ^200), 
((rii + l)ao, (^2 + l)ao)}, respectively (see fig. 1). 

Let us now for a moment take the more common point of view and consider 
the skeleton we have constructed above as a piecewise fiat manifold, and not as 
a discretization of a continuum manifold. The curvature is then concentrated at 
the vertices and is given by J^i^i^i^ — Xi) 0, where 5i is the deficit angle at the 
vertex Xj. The deficit angle 5i is given by 5j = 27r — (pt, where 0t, t = 1, . . . , 6 
are the six dihedral angles that meet at site Xi. From the point of view that a 
skeleton represents an approximation of an underlying smooth continuum manifold, 
such a singular curvature distribution is unnatural. Rather we should approximate 
the curvature in a region Ai around the vertex Xi by a constant Ri, such that 
RiAi = d'^x^J g{x)R{x). This implies that 

Ri = 26i/A, . (2.7) 

There is no unique prescription to assign an area Ai to a vertex Xj G Tq. The 
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easiest way is to follow ref. and use a baricentric subdivision where one divides 
each triangle on T along the side bisections into three equal pieces. The area Ai 
associated with the vertex Xi is then given by 

6 

A = lY.^^t, (2-8) 

t=l 

where Au, t = 1, . . . , 6 are the six triangle areas connected to the vertex i. For a 
given triangle with edge lengths la, Ip and /-y, the triangle area Aap<^ and the dihedral 
angle 0^ which is opposite to the edge are given by 



Ao^p, = + 2(111} + + (2.9) 

With the replacement R{x) — > 26i/Ai for x G it is also straightforward to dis- 
cretize the term which poses a problem on a piecewise flat manifold because 
one cannot square the delta functions. This leads to the Regge form of the action 

(ED. 

S{1) = {XA, + 2k6, + 4(36^ /A,} , (2.11) 

i 

with 2J2i^i = the discrete version of the Gauss-Bonnet theorem. For a given 
torus with a sufficiently smooth metric field, this discretized Regge action approaches 
the continuum action ( p.2|) in the limit that the regularization is removed. It can be 

checked, for instance using the Mathematica package, that Ai = g{xi)al + 0{al) 
and Ri = R{xi) +0(ao), when expressing the edge lengths in eqs. ( |2.7] ) and ( ^.8[ ) in 
terms of the metric, using relations ( p.6|) , and expanding to leading order in oq. 

The above discretization procedure is not invariant under general coordinate 
transformations on T. Let us illustrate this for the cosmological term. Using the 
relations (2^) and (2^) we can express the sum in (2^) in terms of g^^, 



a? 

Ai = al^+ — ^ {-gii d^^g22 - g22dx^gi 



+912(9x^922 + dx^gn) + 2gi2{dx, + dx,)gi2} + 0{a^) , (2.12) 

with the metric and its derivatives evaluated at the point Xi = {rii, 77-2). The relation 
(|2.12 ) shows that the general coordinate or diffeomorphism invariance is broken by 



the 0(ao) terms in the action. This implies that a given set of continuum manifolds 
which are related by general coordinate transformations on T will lead to Regge 
skeletons whose actions differ by 0(ao) effects. 
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There are two sources of discretization effects when translating a continuum 
manifold into a Regge skeleton in the above way: Firstly, of course, discretization 
effects depend on the value of A^. Secondly, the discretization effects depend on 
the choice of coordinates on T. By choosing an optimal coordinate system, one can 
minimize the discretization effects at a given A^. For instance, one could choose a 
gauge fixing condition that minimizes the square of the terms in the action which 
are oc (cf. eq. ( |2.12| )). 

According to the picture sketched above, the edge lengths in a Regge skeleton rep- 
resent both physical degrees of freedom, corresponding to curvature fluctuations, and 
gauge degrees of freedom, corresponding to coordinate reparametrizations. Since in 
general there does not exist a transformation which is an invariance of the Regge 
action ( |2.11| ), the identification of a "gauge transformation" on the edge lengths 
is necessarily ambiguous. For instance, such gauge transformations on / may be 
defined as 

1/2 



I, 



/? + 29;ei 
il + 2d',i2 



1/2 



il/2 



(2.13) 



where d'^^^ = {iu^i+ji — ^u,i)/cio- This transformation is constructed such that after 
inserting relations ( |2.6| ) into ( p.l3| ), one recovers from the leading order terms in ^ 
the infinitesimal continuum gauge transformations, g^j^i, g^u + d^^^ + d^^^. The 
transformations (|2.13|) are associative, i.e. l^^^ = (l^)^ . This is an important 
property since it allows us to devide the configuration space of Regge skeletons into 
equivalence classes (gauge orbits), that consist of skeletons which are related by 
a gauge transform of the form ( p.l3|) . 



2.3 Integration measure 

Next we turn to the problem of defining the measure for the integration over Regge 
skeletons in the path integral. In the continuum path integral, the measure (before 
gauge fixing) is of the form 



x 



(2.14) 



For a 



—1/2 this measure is scale invariant, for a = —1/4 it is gauge invariant 

One could now use the relation (Of) to 



as follows from Fujikawa's method [22 



regularize this measure by expressing it in terms of the edge lengths /, 



g{x 



,3a 



n dgt„y 



X 



d/i(o = |[2A23(o]'"n^^^ 



(2.15) 



f3=l 
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where we have suppressed the site index i. The triangle area A123 is defined in (^.91). 



Such a measure includes the gauge degrees of freedom. In the continuum the gauge 
degrees of freedom are removed by gauge fixing and the full measure would acquire 
an additional gauge fixing and Faddeev- Popov term. For the integration over Regge 
skeletons one is not forced to fix the gauge because the gauge invariance is broken 
by the the discretization procedure and no infinite volume of the diffeomorphism 
group can arise. 

At this point one can follow two strategies: One can ignore the gauge degrees of 
freedom (the ^'s of eqs. ( |2.13| )) and hope that these modes decouple automatically 
from the physical degrees of freedom when approaching the continuum limit. Small 
^ fiuctuation are indeed almost decoupled, when the Regge skeleton is sufficiently 
smooth as has been demonstrated in ref. |]15|, but for a generic skeleton they are 



strongly coupled and it is difficult to beheve that they can ever decouple from the 
physical degrees of freedom in a non-perturbative evaluation of the path integral. 
For a discussion of this scenario of dynamical gauge symmetry restoration in the 
context of a non-perturbative formulation of chiral gauge theories on the lattice, see 



ref. 23 



The second strategy one could follow, is to mimic the continuum approach. This 
implies that the gauge degrees of freedom should be removed from the path integral 
by employing a suitable gauge fixing and Faddeev-Popov procedure. To implement 
this in the Regge path integral, one could translate the continuum gauge fixing term 
and corresponding Faddeev-Popov ghost term into a Regge form, e.g. using the 
correspondence (|2.6|) . In the continuum the BRST symmetry prevents that gauge 
non-invariant terms emerge in the low energy theory; on the lattice this symmetry 
is not present and one expects that gauge non-invariant counter terms have to be 
added in order to recover a BRST invariant low energy theory. Such an approach has 
been advocated some time ago for the lattice regularization of chiral gauge theories 
[Q, which also have the difficulty that the local gauge symmetry is broken by the 
regularization. The technical problem with this approach is that the ghost term in 
the action is very hard to include in a numerical Monte Carlo evaluation of the path 
integral and a non-perturbative computation of the counterterm coefficients may be 
very cumbersome. 

This second option appears to have the best chance of leading to a well-defined 
path integral, since it should work at least in perturbation theory. However, faced 
with the technical difficulties of this strategy, previous studies of the Regge path 
integral have resorted to the first (non-gauge fixing) approach. The integration 
measure was chosen to be of a very simple form 

^/^C(0 = 1171^4, (2.16) 

a.i oa 
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which is scale invariant for C = 0. For a recent investigation of the ( dependence of 



expectation values of local observables see ref. ||25[] . 



3 Regge scaling relations 

Consider now the path integral for Regge skeletons with a fixed area A, 



Z{A, Ni) = J dn^il) F(l) exp 



5(5: A, -A) 



(3.1) 



where we use the abbreviation = iSf /Ai. The total curvature of a skeleton is Anx 
and we leave this term out from the action since we shall study only skeletons with a 
fixed topology. The measure dfi(^{l) is of the simple form ( p.l6| ) and is scale invariant 
for ( = 0. The factor F{1) is equal to zero if the skeleton contains one or more 
triangles for which the triangle inequalities are violated, and is one otherwise. We 
have in ( |3.1| ) included the argument A^i, which is the number of edges (1-simplices) 
in the skeleton, to indicate that we evaluate the path integral for a fixed number of 
edges with fixed connectivities. 

To eliminate the delta function, we first perform a transformation of variables. 



(3.2) 



Since F{1) which imposes the triangle inequalities is unaffected by this change of 
integration variables, eq. ( p.l|) can be written in the form 



Z{A,Ni) = J c//ic(v^/) F{1) exp 



5(3 A -A). (3.3) 



Since the partition function Z{A,Ni) does not depend on s we may integrate both 
sides in ( [3.3[ ) over s with a normalized weight function P{s), i.e. dsP{s) = 1, 



Z{A,Ni) = J dii^{l)F{l) J dss^''^/^P{s)exp 



i s 

' A 
Ai 



exp[-\A] / dfx^{l)F{l): 



A 



JViC/2 



J2i Ai \ Ai 

d/x^/) F{1) exp[-^-^'^(/)] , 



P 



5isJ2A-A) 



exp 



(3.4) 



where the last equality defines the rescaled action S^'^'^^'^ (l) . 

Now it appears as if Z{A, Ni) depends on the arbitrary function P, but this is 
of course not the case. This follows from an infinite set of Ward identities which 
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can be derived by using again the scale transformation ( p.2| ), and observing that 
Ni)/ds'^ = 0, n = 1, 2, 3, . . .. Forn = 1 this gives the identity 




-1 



where the expectation value is with respect to the rescaled action, 
(.) = Z{A,N,r'l . exph^^^-'^(/)] . 



(3.5) 



(3.6) 



We want to investigate if the Regge partition function Z{A, Ni) reproduces the 

oo. Therefore it is useful 



scaling behavior ( |2.4| ) in the continuum limit, i.e. for A^i 
to consider the derivative 



dlogZ{A,Ni 
dA 



(3.7) 



with A/5 and Ni kept fixed. Here we have used the identity (|3.5|) to eliminate a P' / P 
term. 

Next we note that the expectation value Aj J2i Rf) is a function of the ratio 
(3 /A and of A^i. This follows from inspecting the path integral ( p.4| ). If for A^^i oo 
this expectation value behaves such that 



/3 



A ' 2 

then the Regge partition function would scale as 

Z{A, oo) oc exp 



= (5A + (57 + 0(-) 




(3.^ 



(3.9) 



and we may identify A + 6X/j3 with the renormalized cosmological constant and 
+ 2 with the string susceptibility 7str- 

Since the term {l3/A) J2j Aj J2i Rl is the total action in the path integral (|3.4|), 
its expectation value in general will grow proportionally to A^^i. Hence the limit 
Ni ^ oo only exists, if the exponent ( in the measure is chosen such that Ni(/2 
cancels the term oc A'^i in the expectation value of the action. 

Instead of the measure ( |2.16D we could also have used a different measure where 
the exponent ( is given by C = Co + Ci/^i- Now the term Ni(q/2 has to cancel the 
divergence, but the finite part is shifted by the arbitrary constant (i which gives a 
contribution to the string susceptibility. Hence ^7 appears to depend very sensitively 
on the choice of the measure. The only way that 5'y could be independent of Ci is 
if it is compensated by the expectation value of the action: An expansion to lowest 
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order in Ci/^i shows that 7str is shifted by + c) with c the connected part of 
((So.i log^m)(I]j Si /^-Rf/^)); evaluated for (i = 0. It is very unhkely that c is 
exactly equal to —1/2, such that ^7 is independent of Ci- 

This sensitive dependence of ^7 on the measure is seen particularly clearly for 
j3 = 0: Then eq. ( p.7| ) can be integrated easily, giving 

Z{A, Ni) oc A-i+^i^»/2+Ci/2 exp[-AA] . (3.10) 

For the choice Co = the limit A'^i 00 can be taken and the partition function 
is of the same form as in the continuum, but with an arbitrary string susceptibility 
7str = 2 + Ci/2. 

This discussion shows that the real difficulty is the lack of a criterion to fix 
such ambiguities. This is different in the continuum where the power ( is fixed in 
Fujikawa's method [^, by the requirement of gauge invariance. Such a guideline is 
lacking in the Regge model. 

These considerations cast serious doubts on the claim made in ref. |1^ that the 
Regge model with scale invariant measure reproduces the continuum values of the 
string susceptibility. To settle this, we have also performed a numerical simulation 
of this model using the scale invariant measure, i.e. for C = 0, and three different 
topologies with genus h = 0,1 and 2. We shall see that the numerical estimates for 
7str differ substantially from the continuum result for h ^ 1. 

4 Details of the numerical simulation 
4.1 Histogramming 

The area dependence of the Regge partition function can be inferred from eq. (|3.7|) 
by computing {{P / A) J2j J2i n't) = —Pd\ogZ{A,Ni)/dp in the rescaled model 
( |3.4|) . This is identical to evaluating the expectation value of in the original 
model with the delta function (|3.1|). Hence there are two ways to proceed: One can 
simulate either the original model, including the delta function that constraints the 
area to be fixed, or alternatively one can simulate the rescaled model ( |3.4| ). 

The first option poses the problem of including the constraint that the total area 
must be fixed. Such a global constraint is hard to implement efficiently in a Monte 
Carlo simulation (see however ref. ^6]). Alternatively it is possible to include the 
constraint a posteriori. To this extend we first approximate the delta function by a 
narrow block function which is normalized to one. 



(4.1) 
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Then the expectation value of the observable 0{l) is given by, 

(O)^ = ZiA,N,)-'l dfx^il) F(/)0(/) A^,iY.^^-A) exp[-5(/)] 

i 

= (0(/) A,(E. A, - A,(E. A, - A)) , (4.2) 

where the action S{1) is given in eq. ( |2.11| ) with k, = and (•) in the lower line 
is evaluated in the model without the area constraint. In practise this amounts to 
generating an equilibrium ensemble of skeletons with the unconstraint action and 
then computing the average of f3R^ on the subensemble for which the total area is 
within the interval [A — h, A + h] for a selected value of A. We shall refer to this 
method as "the histogramming technique" . 

Most of our numerical results to be presented below, have been obtained with 
this method. For the simulation of the unconstraint action we use a Metropolis 
algorithm. The term in the Regge action extends over several triangles and in 
order to make the code vectorizable we had to use a four color subdivision of the 
lattice. In the model without the area constraint we can derive the Ward identity 



+ ^Jdfi^ F{1) (^Xj2A-(3 J2Rl^ exp[-5(/)] = , 



(4.3) 



2 

Z = J dj^^ F{1) exp[-^(/)] , (4.4) 

which follows from the invariance under rescaling of the links. This shows that 
the unconstraint path integral is divergent for C = ^^^^d (3 = because the ex- 



pectation value of > can only vanish if (^^) diverges. This is perhaps 



seen even more clearly by writing the unconstraint partition function in the form 



Z = dA Z{A,Ni) with Z{A,Ni) given in eq. ( p.l|) . The integration of rela- 
tion (P?7] ) for = yields Z{A,Ni) oc A'^/'^'^ exp[-XA], which shows that the 
unconstraint partition function Z has a non-integrable singularity at A = ii ( < 0. 

Also the rescaled action (|3.4|) could be simulated using a Metropolis algorithm. 
However, the non-local action j3 J2j J2i R^ prevents vectorization of the code. This 
suggests to use a Hybrid Monte Carlo algorithm, which can easily be vectorized and 
is usually found to be more efficient than a Metropolis algorithm. 



4.2 Hybrid Monte Carlo algorithm 

To use the Hybrid Monte Carlo algorithm we first introduce new degrees of freedom 
Pai which are viewed as the canonical momenta of the edge lengths lai with respect 
to the measure dl^i dpai- The path integral ( p.4| ) may then be written in the form 

Z = J Y[{dL^ dp^i) F{1) exp[-H{l,p)] , (4.5) 

a,i 
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H{l,p) = IY.pI + S{1), (4.6) 
+ (|iViC + l) ln[^A,] -(C-l) ElnL^, (4.7) 

\ i / a,i 

where we have introduced the hamiltonian H{l,p) and dropped the constant \A — 
^Ni(\nA from S{1). The various factors in ( p.4|) have been included as logarithms 
in the action, except for F{1) which is a discontinuous function of the edge length /. 

We now want to find trajectories for l{t) and p{t) such that the hamiltonian 
H{l,p) is conserved in the (computer) time t. After choosing the canonical relation 



we can compute the time derivative of p from the requirement that dH/dt 
which leads to 

dpai _ dS{l) 



(4.9) 



So far we have not done anything to include the constraint F{1) on the edge 
lengths. The region in configurations space where the triangle inequalities are vio- 
lated can be avoided by setting the action in this subspace equal to infinity. The 
trajectories then should bounce back from this potential barrier, with the condition 
that the momentum is conserved. One way to achieve this is to flip the sign of the 
momentum component perpendicular to the plane Pa—Pp—p-y = 0, when the triangle 
with edge lengths /q, Ip and just starts to fail the constraint /q, — //j — > 0. A less 
sophisticated procedure is to flip the signs of all momenta, p^ — > — Pa, Pp ~Pf3, 
P7 -Pf 

A new configuration {lai} is computed by integrating the hamiltonian equations 
of motion ([4.8|) and (|4.9|) from ttot+At. These equations are integrated numerically 



by means of the standard leap frog method [27|, which preserves the time reversibility 



of the evolution. To compensate for discretization errors, H{t + At) = H{t) + 
0(5t^), one accepts the new configuration {l(t + At),p{t + At)} with probability 
min{l, exp[— (if(t + At) — H(t))]}. We use a discretization time step 6t = At/n, 
with n ~ 10 and 6t chosen such that the acceptance rate is between 60 and 80%. 

A practical way to implement the constraint F{1) is the following: We check 
after each time step 6t if one or more triangle inequalities are violated. If this is the 
case, we assume that the momenta are constant over the interval St and compute the 
earliest time when such a violation has occurred. We use one of the above mentioned 
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A C 



Figure 2: Triangulation of the sphere with N = 3 (fig. a) and the bi-torus with 
N = 4 (fig. b). The boundaries with the same labels have to be identified. The full 
circles mark the vertices whose coordination number differs from six. 

prescriptions to flip the momenta which correspond to the hnks that are responsible 
for the triangle inequality violation. Afterwards we try to complete the time step 
with the new set of momenta and repeat the above procedure of flipping momenta 
until no triangle inequalities are violated at t + 6t. This procedure, which is easy to 
implement in two dimensions, ensures that the time evolution is still reversible in t 
when the discretized trajectory bounces off the regions in configuration space where 
triangle inequalities are violated. We find that the number of triangle inequality 
violations is typically small and decreases at larger jS. 

From eqs. ( p.9|) and (|4.7| ) it follows that dS{l)/dla contains terms of the form 
dAafB^/dla = la{—la + ^| + ^'^) /^Aa/s-y, which diverge when If^ ^ la + l^ or — > 
la + 1/3- This implies that discretization errors may blow up near the excluded 
regions in configuration space, which implies that the time step 6t has to be taken 
very small in order to keep the acceptance rate sufficiently large. The reduction of 
6t makes the algorithm relatively slow. We could improve on this by generalizing 
the triangle inequalities, F{1) — >■ F^{1), such that the minimum area that is allowed, 
is increased from zero to some small non-zero value, i.e. by using the inequalities 
—la + 1/3 + 1^ > e > instead of —la + lf3 + l^ > 0. With this modification the Hybrid 
Monte Carlo algorithm becomes competitive in speed with the Metropolis algorithm 
and is more efficient for larger values of f3. 
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4.3 Topologies with genus zero, one and two 

Let us close this technical section with a description of the various topologies we have 
used in our simulations. The simplest skeleton has the topology of a torus, cf. fig. [|. 
We use the same simple triangulation as in ref. |T^, which is obtained by taking 
a regular lattice with periodic boundary conditions in spatial and time directions 
and adding diagonal links, such that each site x = {nia,n2a) is connected by six 
links, as discussed in sect. 2.2. On a torus with A^o = vertices (0-simplices), the 
number or triangles (2-simplices) is N2 = 2Nq and the number of links (1-simplices) 
is Ni = 3Nq. The Euler characteristic is x = -^2 ~ -^1 + -^0 = 0. 

Our second topology is that of a sphere, which has the Euler characteristic 



X = 2. Here we do not follow the construction of ref. |]T6|, in which a cylinder is 
made into a sphere by pinching each of the two boundaries to a point. These two 
vertices have an exceedingly large coordination number of order -y/iVo which has 



been claimed in ref. [0 to affect the scaling of the partition function. To avoid 
large coordination numbers, we start from a three dimensional cube and use the two 
dimensional lattice consisting of its boundary. Each of its six faces is triangulated 
in the same regular way as the torus. In fig. ^ we have cut the cube open in order 
to show the triangulation of its front and back side. All vertices have coordination 
number six except for six of the eight corner points whose coordination number is 
four. If the number of sites on the square faces of the cube is iV^, the total number of 
vertices on the sphere is Nq = 6(A^ — 1)^ + 2, the number of triangles is N2 = 2No — 4: 
and the number of links is A^^i = 3No — 6. 

To illustrate how typical surfaces look like, we have embedded the two dimen- 
sional surfaces of the sphere in a three dimensional space. For typical equilibrium 
configuration at a small and large value of the coupling (3 /A we have determined 
the embedding space coordinates X/, Xf and Xf associated with vertex i using the 
relations /^^ = I]^=i(Xf — XjY, with i and j denoting the two neighboring vertices 
connected by the edge lai- Fig. |^ exhibits two examples for such surfaces, one ob- 
tained at j3/A = 0.5 (left graph) and the other at P/A = 0.0005 (right graph). It is 
clearly seen that the surface at P/A = 0.5 is smoother and appears to be much closer 
to a sphere, whereas the surface at f3/A = 0.0005 is more irregular and crumbled. 
Notice that the underlying features of the cube have completely disappeared and 
parts of the surface which contain a vertex with coordination number four do not 
to carry an exceptional amount of curvature. 

Finally we have also constructed skeletons with x = ~2, which have two handles. 
These skeletons are obtained by gluing together two tori of size N"^, as illustrated 
in fig. 03. In both tori we cut out a square window of size A^/2 and glue them 
together by identifying the links that are on the boundaries of the two windows. 
The four corner points in the window frame have a coordination number larger 
than six: two have coordination number eight and two have ten. All other vertices 
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Figure 3: Equilibrium surfaces with the topology of a sphere, embedded in three 
dimensional fiat space. The couplings are (3/ A = 0.5 (left graph) and 13/ A = 
0.0005 (right graph). The number of triangles is N2 = 192. In the left graph one 
vertex with coordination number four is visible. 

have coordination number six. The total number of vertices on such a bi-torus is 
Nq = 3N'^/2 — 2 {N must be even), the number or triangles is N2 = 2Nq + 4 and 
the number of links is A^i = 3Nq + 6. 

5 Numerical results 

In order to compute the string susceptibility and renormalized cosmological constant 
we have to compute the A/f3 dependence of the average action in the constraint 
area model for N2 — >■ 00. This term should behave as {P/A){J2jAj J^iRf) = 
(A/P) SX + + 0{j3/A), which leads to the string susceptibility 7str = ^7 + 2, cf. 
eqs. ( |3.8| ) and ( |3.9| ) above. 

Most of our data have been generated with the unconstraint action. We used the 
histogramming method to compute the average action (3{J2i RI)a with the constraint 
expectation value {•) defined in eq. (|4.2| ). For a fixed value of (3 the distribution 
of triangle areas in an ensemble is peaked around an average triangle area a = 
Ai) / N2, with A^2 the number of triangles in the skeleton. For increasing N2 the 
peak of the distribution becomes more narrow, but remains roughly at the same 
value of a. In our numerical simulations we have kept (3 fixed when increasing N2. 
Eq. ([4.3|) shows that it then would require tuning (3 (x N2 to keep the bin with a 
fixed value of Y^i^il (3 sufficiently populated when increasing the skeleton size. 
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It requires several steps to estimate a value for the shift of the string susceptibility 
S^: First we generate a large ensemble of 10^ — 3 x 10^ configurations for each 
lattice size and store for each configuration the values of J2i and J2i -Rf • Then 
we use the histogramming technique to compute the expectation value Sji2{/3/a) = 
l3{Y^i R'i/N2)j^ for several values of the ratio P/a. In fig. | we have plotted Sr2[(3 /a) 
as a function oi (3/a for skeletons with topologies of a bi-torus (fig. §a), a sphere 
(fig. ^) and a torus (fig. ^c). The circles indicate the centers of the bins. On 
the larger lattices the width of the distribution is substantially larger than the (5 /a 
interval shown in fig. ^ and most of the data points are outside the frames. Each 
figure contains the results for a number of different lattice sizes. These distributions 
overlap over the full range of values of N2 which makes it easy to read of the N2 
dependence. As noted above, this would not be the case if we would plot Sr2 as a 
function of 13/ Ai. The solid lines in figs. b, and c have been obtained by a 
polynomial fit which has been used to interpolate these results to obtain 6*^2 (/3/a) 
over a range of (3/a values. The fits on the smaller lattices include also those data 
which are outside the frames. 

In figs. b and c we have displayed the interpolated values for S'ij2(/5/a) as a 
function of I/N2 for the three different topologies and several fixed values of (3 /a. 
It is seen that for the larger values of N2 shown in the figure, the data fall nicely 
on straight lines. For smaller N2 we find deviations from the straight line behavior. 
This shows that for sufficiently large N2 Sji2[f3/a) behaves as 

SR2{f3/a) = Co{(3/a) + c,{(3/a)/N2 + ■■■ . (5.1) 

The coefficients cq and ci are found using a fit (straight lines in figs. b and 
c). The value of was in all cases smaller than two. 

Next we have to convert the results for cq and Ci into an estimate for the string 
susceptibility and the renormalized cosmological constant. In sect. 3 we have shown 
that, 

{f3E^ R^)a = N2 Sr2{N2(3/A) = N2 co(iV2/?M) + ci{N2(3/A) + • • • (5.2) 
should behave as, cf. eq. ( |3.8| ), 

{PT.^ R')a = - I + 57 + OiP/A) . (5.3) 

This is the case if the coefficients Cq and Ci can be expanded as Co,i(a;) = Cg^l + 
Cai/^ + 0{l/x'^). Combining all results then leads to the relations 



(5.4) 
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In order to cancel the term oc A^^i in the right hand side of eq. ( |3.8| ), we should 
choose the measure such that 

C + AcP/3 = 0. (5.5) 

Here we have used the relation N1/N2 = 3/2 which is valid for any two dimensional 
skeleton. 

In order to make contact with ref . ||TB| , we have used the scale invariant measure, 



( = 0. Then relation ( p.5|) is not satisfied, but we shall nevertheless define a string 
susceptibility from the finite part of {J^iRDa' 

i.e. 7str = cf^ + 2 = ci(oo) + 2. We 
shall show below that cf^ in fact appears not to change much when tuning ( such 



that dJ) is fulfilled. 

To determine this susceptibility from our numerical data we have to extrapolate 
the Ci(/?/a) data to P/a = 00. The results for Ci(/?/a) are plotted in fig. ^ as a 
function oi (3/a for the bi-torus (fig. a) and the sphere (fig. b). The two plots show 
that it is difficult to extrapolate reliably to large values of [3 /a. For the bi-torus 
we find from our data a lower limit ^7 = c^^^^S.S, provided that ci{(3/a) keeps 
monotonously increasing for (3 /a > 15. Using the heuristic fit ansatz 

Ci{x)=A + B/{x + C) (5.6) 

we find plausible values for in the range 4.0 ~ ~ 5.0. The curves obtained 
by fitting the ci(/3/a) data to the ansatz ( pl6|) are represented in fig. |]by the solid 
lines. 

For the sphere we have computed Ci (/5 / a) only for a few P/a values. Fig. ^jb shows 
that the /3/a dependence of the coefficient Ci{j3/a) is much weaker than for the bi- 
torus, but also in this case ci{(3/a) appears to increase monotonously with {3 /a. The 
data which are represented by the full circles have been obtained with the fixed area 
Hybrid Monte Carlo code. We have used the generalized triangle inequalities which 
put a non-zero lower limit on the minimum triangle area (e = 0.01), as discussed 
in sect. 4.2. As can be seen in fig. such a small value of e appears not to lead 
to significant deviations compared with the Metropolis results. From this figure we 
can read of a lower bound for cf' which is approximately equal to 3.5. 

These results were obtained using the same scale invariant measure as in ref. 
1^ . It is clear that the resulting values for 7sti. = 2 + are in strong disagreement 



with the continuum predictions for both topologies: For the bi-torus we find 7str~5.5 
whereas the continuum value is 7str = 4.5 and for the sphere we find 7str~5.5 whereas 
the continuum value is 7str = —0.5. The continuum values of 7str — 2 are represented 
in fig. ^ and b by the horizontal dashed lines. For the torus we find that Ci {{3 /a) ^ 0, 
independent of (3 /a and hence also ^7 = cf'^ ~ 0, as in the continuum. This means 
that the partition function obeys the scaling relation Z{A) oc A~^, which, as we 
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pointed out in sect. 2.1, also holds for the classical model. Therefore it gives no 
information on the contribution of quantum effects, and finding the correct scaling 
of the Regge partition function for this case is of little significance. 



The results for the bi-torus and sphere show that the claim of ref. [|16[ , that the 
Regge model with a scale invariant measure reproduces the continuum values of the 
string susceptibility, is not justified. However, we have argued above that the value 
of ( in the measure has to be chosen such that the term oc Ni in the expectation 
value of the action is canceled, cf. eq. ( pl5|) . This is not the case when using the 
scale invariant measure: e.g. we find for the bi-torus co{(3 /a) = 0.1072(4), 0.1627(7), 
0.2096(2), 0.2225(2), 0.2316(2), 0.2355(6) for p/a = 0.192, 0.61, 2.4, 4.5, 8.6, 12.6 
and clearly ( |5.5| ) is not fulfilled since C = in our case. 

To check if the numerical results change when adjusting ( such that cq + 3^/4 
becomes closer to zero, we have performed a simulation at /3 = 0.1 and ( = —0.2 
for the case of the bi-torus. For (3/a = 1.7 we find cq + 3^/4 ~ 0.0904(4) which is 
indeed almost by a factor two smaller than the value of cq + 3^/4 obtained for C = 
at the same value of [3/a. The result for the coefficient ci is represented in fig. ^ 
by the triangle. Even though the error bar is rather large, it occurs that the result 
falls on top of the fit curve to the C = results which indicates that the coefficient 
ci(/3/a) does not depend on the measure parameter C^. 



6 Summary and final remarks 

In this paper we have investigated the Regge method of regularizing the euclidean 
path integral for quantum gravity, using the two dimensional model as a test case. 
We have shown that a Regge skeleton can be viewed as a discretization of a smooth 
continuum manifold, which is defined by a metric field g^y. The Regge action, 
including terms, approximates the corresponding continuum action, up to dis- 
cretization errors. These discretization errors break diffeomorphism invariance, as 



illustrated in eq. ( |2.12| ) for the cosmological term in the action. 

This suggests to separate the Regge degrees of freedom, the edge lengths of 
the skeleton, into equivalence classes, as in the continuum. We have given such a 
possible prescription in which skeletons are equivalent if they are related by a "gauge 
transformation" / l^, cf. eq. ( |2.13| ). In a Gaussian approximation of the action 



around fiat space, these transformations are zero modes, and correspond to the 
continuum gauge transformations [|l^. The presense of these modes suggests that 
the measure in the Regge path integral should include a gauge fixing and Faddeev- 
Popov term, as in the continuum. However, it is not known how to implement this in 
practical simulations, and so far only simple variations on a scale invariant measure 
have been used. 

Our numerical simulations have been done with such a simple measure. One can 
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derive scaling relations for the Regge path integral of two dimensional gravity, 
which relates d\ogZ{A)/dA to the expectation value of the curvature square term, 



evaluated at a fixed value A of the total area of the two dimensional surface |jT6[ . We 
use two Monte Carlo methods to compute the constraint expectation values. In the 
first of these methods we simulate the unconstraint path integral and implement the 
fixed area constraint a posteriori by using a histogramming technique. The second 
method is a Hybrid Monte Carlo algorithm with which such a constraint model 
can be simulated directly. This algorithm is very attractive, because it can easily 
be vectorized and in comparison to the histogramming method it is more efficient 
because it allows us to carry out a simulation at a chosen value of (3 /a. 

From the exact scaling relations for the Regge partition function we infer that 
the string susceptibility is determined from the part of P{R'^)^ + Ni(/2 that is finite 
after removing the regularization, i.e. for A^^i oo. The first term is the expectation 
value of the action at a fixed value of the area A and ( is the measure parameter. 
The identification of this finite part is however ambiguous and depends sensitively 
on the specific form of the measure used, as can be seen by shifting the measure 
parameter ( ( + (i/Ni. This is particularly clear for /5 = 0, where we find the 
exact result 7str = Ci/2- 

These considerations cast serious doubts on the claim made in ref. that the 
Regge model with scale invariant measure reproduces the continuum values of the 
string susceptibility. To settle this, we have also performed a numerical simulation 
of this model: as improvements over the earlier work, we have done the simulation at 
a fixed value of the total area, we have used a spherical topology with a more regular 
skeleton which only has coordination numbers four and six, and we also included a 
topology with genus h = 2, the bi-torus. Otherwise we use the same scale invariant 
measure and R^ term as in ref. |T^. We have determined the coefficient Ci{(3/a) of 
the I/N2 term in the average action to a high accuracy for all three topologies and 
for several values of the coupling constant [3 /a. 



In ref. the coefficient Ci was wrongly identified with the string susceptibility. 



It turns out that Ci depends strongly on j3/a for the case of the bi-torus. For the 
sphere the (5 /a dependence is much weaker (see fig. D and for the torus we find no 
(3 /a dependence. We have argued that the determination of the string susceptibility 
requires that P/a ^ 00 and that ( is tuned such that eq. ( p.5[ ) is satisfied. In our 
simulations we have followed ref. |]16| and chosen the scale invariant measure C = 0; 
for which eq. ( |5.5D is not fulfilled. Simulations at a negative value of ( at which this 
equality is almost satisfied indicate that ci is not affected by a small change in (. 
Assuming that we do not include a term (i/Ni in the measure, we can then from 
the Ci(/5/a) values shown in fig. ^ infer lower bounds on 7str- Both for the sphere 
and the bi-torus we find 7str~5.5 which is substantially larger than the continuum 
values given in eq. (p75|). Only for the torus the string susceptibility reproduces the 
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continuum value, but here the naive scaling holds, which does not involve quantum 
fluctuations. 

A likely explanation for this failure of the Regge approach is the contribution of 
the gauge degrees of freedom, the ^'s of eq. ( |2.13D . These modes do not decouple 
form the action, because the Regge action is not gauge invariant (this is only approx- 
imately the case for small transformations of a smooth skeleton |]T^). The gauge 
modes have also not been removed by including a gauge fixing and Faddeev-Popov 
term in the integration measure. If this is indeed the source of the failure of the 
Regge approach in two dimensional gravity, one would expect similar problems also 
in four dimensions. 

To test this hypothesis and to improve the Regge approach, one could try to 
discretize the continuum model, including these gauge fixing and Faddeev-Popov 
terms. Such an approach faces many complications: The main obstacle is that one 
would have to find a practical method to simulate an action that includes such a 
gauge fixing term and Faddeev-Popov ghosts. Another serious complication arises 
from the breaking of the gauge (diffeomorphism) invariance which will produce gauge 
variant terms under renormalization. These terms have to be canceled by adding 
suitable counterterms. Furthermore, such gauge fixing and Faddeev-Popov ghost 
action expressed in terms of the Regge edge lengths, will be very cumbersome. Then 
it might be easier to abandon the Regge approach and use a more straightforward 
regularization of the path integral by discretizing the metric field g^y on a hypercubic 
lattice. 
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Figure 4: Histograms for the bi-torus (fig. a), the sphere (fig. b) and the torus 
(fig. c). We have plotted SR2[j3/a) = P {J2i Ri / as a function of (3/a for several 
lattices (bi-torus: from top to bottom N2 — 108, . . . , 1728; sphere: from top to bottom 
N2 = 192, ...,4332; torus: from bottom to top N2 = 18, . . . , 4608;. The circles 
represent the data points. We have dropped the error bars when they are smaller 
than the symbol size. 
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Figure 5: SR2[j3/a) = (3{Y.i^/N2) a a function ofl/N2 for the bi-torus (a), the 
sphere (h) and the torus (c). We have not included in these graphs the Sb?{P/o) 
results on small volumes since they exhibit deviations from the straight line behavior. 
The symbols and lines in figs, a-c correspond to various values of (3 /a: (a) □ : (3 /a = 
0.61, o: (3/a = 2.4, A: f3/a = 4.5; (b) U:(3/a = 0.61, •■.(3/a = 2.5, A : /?/a = 4.5,- 
(c) A : (3 /a = 2.4. The full circles in fig. b have been generated with the Hybrid 
Monte Carlo algorithm. The straight lines have been obtained by fitting the numerical 
data to the ansatz ( \5. ]\ ). 
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Figure 6: The coefficient ci{[3/a) as a function of 13 /a for the hi-torus (fig. a) and 
the sphere (fig. b). Most of the numerical results have been obtained with the his- 
togramming method (open circles and triangle). The results represented in fig. b by 
the full circles have been obtained with the fixed area Hybrid Monte Carlo code. The 
full lines are fits with the heuristic ansatz ( \5. dj ). The fit in fig. a does not include 
the five points at (3 /a ~ 0.20. The dashed lines represent the 7str — 2 values in the 
continuum. The triangle in fig. a has been obtained at f3/a = 1.7, ( = —0.2 which 
is closer to the point where l\5.di ) is fulfilled. 
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